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The current simulation for the Capsule Parachute Assembly System (CPAS) lacks fidelity 
in representing added mass for the 116 ft Do ringsail Main parachute. The availability of 3-D 
models of inflating Main canopies allowed for better estimation the enclosed air volume as a 
function of time. This was combined with trajectory state information to estimate the 
components making up measured axial loads. A proof-of-concept for an alternate simulation 
algorithm was developed based on enclosed volume as the primary independent variable 
rather than drag area growth. Databases of volume growth and parachute drag area vs. 
volume were developed for several flight tests. Other state information was read directly from 
test data, rather than numerically propagated. The resulting simulated peak loads were close 
in timing and magnitude to the measured loads data. However, results are very sensitive to 
data curve fitting and may not be suitable for Monte Carlo simulations. It was assumed that 
apparent mass was either negligible or a small fraction of enclosed mass, with little difference 
in results. 


Nomenclature 
BEA = Best Estimate Atmosphere 
BET = Best Estimate Trajectory 
(CrS)p = Parachute drag area 
(CpS)v = Effective drag area of payload or test vehicle 
CDT = Cluster Development Test (series) 
CFD = Computational Fluid Dynamics 
CPAS = Capsule Parachute Assembly System 
Do = Nominal parachute diameter based on constructed area, D, = /4-S, / 1 
DSS = Decelerator System Simulation 
EDU = Engineering Development Unit 
F = Tension force in a parachute riser 
FAST = Flight Analysis and Simulation Tool 
FSI = Fluid Structure Interaction 
g = Acceleration of Earth Gravity 
G = Load factor 
y = Flight path angle 
GPS = Global Positioning System 
h = Height 
HD = High Definition (camera) 
Kae = Ratio of apparent mass to enclosed mass such that ma = me(1 + Kae) 
L, = Reefing line length 
L; = Suspension line length 
ma = Added mass, ma = Me + Ma 
Ma = Apparent mass 
Me = Enclosed mass 
Mp = Mass of parachute and suspension equipment 
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Mass of test vehicle, not including canopies and suspension equipment 


= Over-Inflation Control Line 


Permanent Reefing Line 
Dynamic pressure, q = ; -p- var 


Freestream dynamic pressure 


Humidity-corrected atmospheric density 

Ramp Clear (usually chosen as start of test) 

Serial Number 

Parachute Canopy open reference area based on constructed shape 
Projected frontal canopy area 

Synchronized Position Attitude & Navigation 

Elapsed time 

Tilt angle of the canopy axis relative to the chase aircraft camera 
Inertial velocity 

Total airspeed relative to air mass 

Dry weight of parachute 


Suspended weight of test vehicle, not including canopies and suspension equipment 


I. Introduction 


MONG the largest components of the inflation loads of the Capsule Parachute Assembly System (CPAS) 116 ft 
D, ringsail Main parachutes are the so-called “added mass” effects. CPAS inflation loads are predicted using the 
Flight Analysis and Simulation Tool (FAST), which is anchored by matching peak loads via flight test reconstructions. 
FAST reconstructions of the Main parachutes generally match the peak load magnitude, but the simulated peak load 
timing usually occurs much earlier than the actual data, especially during the disreef to full open. An example for 
Cluster Development Test (CDT)-3-7 is shown in Figure 1. 
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Figure 1. Typical FAST simulation (dashed lines) of disreef to full open with peak loads earlier than test data 


(solid lines). 
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This discrepancy implies a lack of fidelity in the added mass model. Therefore, an effort was undertaken to better 
understand Main parachute added mass. Added mass can be divided into the “enclosed mass” of the air within each 
canopy plus the external “apparent mass” due to viscosity. The enclosed mass is the product of the enclosed air volume 
and ambient air density. 

Fortunately, a wealth of CPAS Main deployment visual and digital information was available to estimate enclosed 
mass. On every flight starting with the Engineering Development Unit (EDU) test series, at least one chase helicopter 
was staged at an altitude to allow for nearly orthogonal side views of the inflating Main parachutes. The payloads 
incorporated upward-looking video cameras to observe the inflating shapes and cluster behavior. The avionics system 
on the payload provided highly accurate system state data,' and the Main risers were instrumented for loads.” 

This paper focuses on identifying the drag and momentum terms that make up the total inflation load to develop 
an alternate inflation model. This model leverages the high-quality data available to CPAS by re-parameterizing the 
system state in terms of enclosed volume. A companion paper? explains how still photos and video were used to 
estimate canopy time-varying volume by constructing 3-D models. An example of these models during the CDT-3-7 
disreef to full open is shown in Figure 2. 


Figure 2. 3-D models of inflating canopies on CDT-3-7 compared to original chase helicopter photo and 
upward-looking HD video still. 


Because calculating added mass requires knowledge of both the amount and rate of change of air mass in the 
canopy, it was important to have a series of photos taken in quick succession. Early in the test program, the 
photographers took fewer photographs during the Main parachute phase. The photographers were later encouraged to 
take photographs at a higher rate during key deployment and disreef events, especially as the cameras improved and 
their storage capacity expanded. A method was developed to synchronize chase photo timing, referenced to an absolute 
scale using the image metadata, and a comprehensive video timeline of known events (such as disreefs and 
touchdown).* A summary of chase photo coverage and the status of volume analysis for Cluster Development Tests 
to date is presented in Table 1. The final three tests in the EDU series incorporated Over-Inflation Control Lines 
(OICL) or Permanent Reefing Lines (PRL) in an effort to control pendulum motion. Because the final design does not 
include such skirt restrictions, those tests are of lower priority for study.° 
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Table 1. Summary of CPAS Main Volume Analysis 


: Main Volume Analysis Status 
Test Vehicle Bay B Bay C Bay E Notes 


Gen II with Added Porosity 
aaa eaEs| 
Sparse early; Chase helo above Mains 
Bay E obscured in 1“ stage 
aaaaicaibi 
ee eee 


Certain features of the 3-D models were validated with previous photogrammetric analysis performed by the JSC- 
KX Image Science & Analysis Group using upward-looking High Definition (HD) videos. That analysis used 
automated tracking to characterize individual canopy skirt perimeter, maximum projected area (S,), and cluster fly- 
out angle (0).7 

A previous study of a single Main inflation was conducted by Airborne Systems using a ground-based camera to 
estimate the canopy volume by approximating an ellipsoid shape.* That study indicated that the added mass effects 
were negligible in the first stage, were more significant in the second stage, and were most significant during the 
disreef to full open, which also matches the current observations. Since the Airborne study, the CPAS Main canopy 
geometry has since changed with the addition of a ring gap and removal of several panels to improve cluster 
performance.’ The current effort also differs from that study by characterizing cluster effects, which often cause the 
canopies to take irregular shapes due to mutual aerodynamic interference. 

The equations for added mass and the simulation methodology is discussed in Section II. Selected flight test results 
are presented in Section II. 


II. Added Mass Equations Using Volume Method 


The current parachute inflation models used by NASA can trace their heritage to a 6-DOF code used by the Viking 
program called UD233.!° The inflation model, including the added mass algorithm, was later utilized by the 
Decelerator System Simulation (DSS).'! During the CPAS program, DSS was modified to better align with the 
DCLDYN (Decelerator Dynamics) simulation, then used by Airborne Systems. !? The DSS inflation parameters were 
modified to allow for common inputs anchored to flight test reconstructions. An added mass filter was added to DSS 
which emulated DCLDYN, but did not incorporate any of its code directly. Unfortunately, there are significant gaps 
in the known heritage of the legacy added mass model. An examination of the DSS source code implies that the 
algorithm was changed over time to match data from available test programs, using geometric approximations and 
undocumented scale factors. When CPAS transitioned to FAST, the algorithms were migrated essentially unchanged. 

The history of added mass has been presented by many authors over the years. Much of this history was recently 
summarized in Ref. 13. The equations for all inflating parachutes are derived from Newton’s second law of motion, 
the general form of which is given in Eq. (1). A simplified free body diagram of a decelerating parachute system is 
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shown in Figure 3. It is assumed that the parachute deceleration vector, the freestream air (Vair) and inertial velocity 
(v) vectors are co-aligned. The elastic spring force is neglected in this study but is included in FAST. 


Parachute Drag: 


F, =4-(Cp), 
Apparent mass, m, 
Enclosed mass, m, 
W,-siny Wp =mp-g 


Riser tension 
force, F 


Vehicle Drag: 


-Flight path | qx (CpS)y 
angle, -v 
Wy-siny 
Var Wy =Mv's 


Inertial velocity, v 


Figure 3. Deceleration with added mass (neglecting spring force). 
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The rate of change of momentum can be written in terms of variable mass multiplied by the system velocity, and 
fixed masses multiplied by accelerations. The relevant terms from the free body diagram of the isolated parachute are 
listed in Eq. (2). 


F+W, sin(—7) ~ Drag=<(m)-v +2 (v)-m (2) 


The items are expanded and re-arranged to yield the components of riser tension (F) in Eq. (3). Added mass (ma) 
is defined in this formulation as the sum of the enclosed air mass (me) and apparent mass (m,). The dry mass of the 
parachute (m,) is multiplied by gravity (g) to produce the weight of the parachute (W,). The CPAS Main canopy 
weight is currently about 250 pounds, including suspension equipment, so the dry mass terms are much smaller than 
the other inflation forces. 


pesg soe ~ 


at V+, +m, + ti, )>-——+ W, siny 


(3) 


A generic block diagram of a typical algorithm that computes riser tension from inflation loads, used by 
simulations such as DSS or FAST, is shown in Figure 4. Inflation parameters such as fill distance constant, exponential 
shape term, and drag area are reconstructed from flight tests. These parameters are inputs to the algorithm and are 
often dispersed, as in Ref. 14. They describe parachute drag area growth with time. Simulations used by CPAS assume 
that drag area monotonically increases with time, although data often show drag reduction in “lagger” canopies due 
to cluster interference. Drag area growth is then related to the growth of canopy volume through some sort of 
geometrical assumption. Drag area is multiplied by dynamic pressure to produce a drag force. The growing air mass 
in the canopy is used to estimate added mass and added mass rate terms, which are then combined with the other 
forces to produce the total parachute riser tension. 
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Figure 4. Typical algorithm for simulating riser tension (e.g. DSS and FAST). 


The formulation of the proposed algorithm, shown in Figure 5, differs in that the enclosed volume state, rather 
than the drag area, is the primary independent variable that is matched to time varying flight test data. The examples 
in this study fit volume data with third order polynomials, which can allow brief volume reduction in a cluster. The 
canopy volume is used to compute drag area using a look-up table that is also derived from flight test data. The canopy 
volume and drag area states are then used to compute riser tension similarly to the typical algorithm. 
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Figure 5. Alternative riser tension algorithm (volume method) with differences highlighted. 


The rate of change of velocity at the vehicle, d(v)/dt, can be calculated using three different methods. Note that 
this term should not include the gravity acceleration vector (typically measured in most accelerometers), because it is 
assumed that enclosed air is neutrally buoyant. CPAS currently instruments its test vehicles with the NovAtel SPAN- 
SE (Synchronized Position Attitude & Navigation),!° which combines GPS readings with an accelerometer housed in 
an IMU'* into an integrated state solution via a Kalman filter. The first method is to directly obtain the inertial change 
of velocity as output from the post processing software, Inertial Explorer.'’ This output was confirmed with two 
independent alternate methods, which may be used if the direct output is not available. One method is to numerically 
differentiate the velocity state in all three directions. The other is to transform the raw acceleration readings from the 
body axis to an inertial frame using Euler angles, and then remove the gravity term from the vertical direction. 

Construction of the drag area look-up table for the proposed volume method involves isolating drag from the flight 
test data through process of elimination. Two other methods were considered but were ultimately discarded. One 
option would be to compute drag through Computation Fluid Dynamics (CFD). The “Bi-model” of Potvin and 
McQuilling computes canopy drag using a series of static CFD solutions for slowly-inflating canopies.'® The Bi-model 
is so named because it is part-way between traditional low-fidelity modeling and Fluid Structure Interaction (FSD. 
Although initially considered for this current study, the Bi-model takes advantage of geometrical assumptions, such 
as radial symmetry and pressurized shape that are not well suited for the clusters of CPAS Mains, which deform due 
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to cluster effects. The other option would be to model the parachutes with FSI. The Team for Advanced Flow 
Simulation and Modeling’? at Rice University, led by Prof. Tayfun Tezduyar, has been successful in modeling many 
features of CPAS-like ringsail parachutes using FSI. These include cluster effects”° and reefed stages.”! However, the 
high computational cost for FSI was prohibitive for this study. 

While computing enclosed mass is fairly straightforward using canopy volume and density data, apparent mass is 
much more difficult to quantify. Potential flow theory provided apparent mass formulations for certain shapes such as 
spheres and ellipsoids. Heinrich used wind tunnel test data and geometric assumptions to postulate that a linearly- 
growing apparent mass term would eventually be equal to 3/8 of enclosed mass by full open, as quoted in Ref. 22. He 
and Noreen noted that even if apparent mass was initially small, the time derivative could be significant.7? Experiments 
by Yavuz and Cockrell indicated that the apparent mass coefficient is more significant than potential flow theory 
would suggest, and varies with angle of attack and an acceleration modulus.™* 

Because this study carefully tracks the growth of enclosed mass, it was assumed that apparent mass could be 
directly related to enclosed mass via a scale factor, Kae. However, initial reconstructions showed an acceptable match 
of peak inflation loads by setting this scale factor to zero. It is likely that any actual apparent mass is therefore being 
erroneously accounted for in either the drag area or enclosed mass terms. Better identifying drag area through CFD 
might reduce this uncertainty. The actual apparent mass formulation may eventually be determined by simultaneously 
optimizing the reconstructions of several tests over a wide range of test conditions. The reconstructions presented 
below initially assume Kae is zero, and each conclude with a sensitivity analysis. 


III. Selected Flight Data Results 


Test analysis begins by validating geometrical inputs to gain confidence in the volume data. During reefed stages, 
the maximum perimeter of the skirt is constrained to be no larger than the length of the reefing line, L;. Each of the 
flights examined here had full open skirt perimeter and projected area photogrammetric data analysis performed by 
the KX Image Science & Analysis Group. As will be shown, these data provide a useful proxy for volume when chase 
photos are scarce. The inflation load components are built up term by term and compared to the measured total riser 
tension. 

A simulation for each canopy is performed at each analysis step. However, rather than constructing a stand-alone 
3-DOF or 6-DOF simulation, the current simulation merely assumes a perfect match of the state by interpolating the 
known Best Estimate Trajectory (BET) data with time. Eventually, the volume method may be incorporated into FAST 
as an alternative inflation model. 


A. CDT-3-2 


The second EDU test, CDT-3-2, descended under two Mains. 3-D surfaces were generated for both canopies and 
are compared to the original images for all three stages in Figure 6. 
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Figure 6. CDT-3-2 selected 3-D canopy models at all three inflation stages. 


The estimated skirt perimeter matches up with the reefing line lengths and KX analysis in Figure 7. 
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This agreement should instill confidence in the computed enclosed volume data plotted in Figure 8. 
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Figure 8. CDT-3-2 enclosed volume data and polynomial curve fits. 


The enclosed volume for each canopy was fit with a third order polynomial for each stage. Main S/N 5 shows 
some reduction in volume at the end of first and second stages. This is consistent with the photos which show that 
Main S/N 6 is the dominant canopy, reducing the airflow into its neighbor. The last data to be fit were immediately 
after inflation, when each canopy finished over-expanding. This feature differs from the current CPAS inflation 
assumptions where each canopy inflates directly to its average full open drag area with no oscillation. The final 
oscillation may be suited to hand off to the CPAS time varying rate of descent model.”° 

It was observed that the enclosed volume makes a linear trend with projected area, as shown in Figure 9. This was 
an unexpected result based on the canopy geometric assumptions made by most inflation codes. For example, the 
volume of a hemisphere is proportional to projected area to the 3/2 power. The shape for a conical frustum/hemisphere 
model (such as that used by Heinrich) would make an S-curve. The slope for each linear fit is shown in the legend. 
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Figure 9. CDT-3-2 linear fit of enclosed volume vs. projected area. 


The assumption of a direct relationship between enclosed volume and projected area allowed for the transformation 
of the volume polynomial curves to projected area curves, which are compared with the original data in Figure 10. 
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The original volume data were numerically differentiated to compute enclosed volume rate through finite 
differencing. Results are plotted in Figure 11. This process is always susceptible to error when there is noise in the 
data. The simulation differentiated each third order volume polynomial, so the results are smoother. Due to the 
approximate nature of the volume curve fits, the simulated volume rates do not match the peaks in the data. This is 
acceptable because those peaks may be exaggerated from the numerical differentiation. Any negative volume rate was 
set to zero, so as to ensure each canopy contributes to system deceleration but not acceleration. 
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Figure 11. CDT-3-2 enclosed volume rate. 


The enclosed mass was computed by multiplying enclosed volume by ambient atmospheric density from the day 
of flight Best Estimate Atmosphere (BEA). Density was stored as a function of altitude, which is interpolated from 
the vehicle state data. The enclosed air mass is shown in Figure 12 and the enclosed air mass rate is shown in Figure 
13. 
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The enclosed mass rate is then multiplied by the vehicle velocity (interpolated from the SPAN-SE state data) to 
compute the added mass rate-based term, plotted in Figure 14. As mentioned previously, although the code carries the 
apparent mass parameter (m,), it is currently assumed to be zero. 
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Figure 14. CDT-3-2 mass rate term. 
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The mass values are multiplied by the derivative of inertial velocity (interpolated from the SPAN-SE state data) 
to compute the added mass velocity rate term, which is plotted in Figure 15. 
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Finally, the parachute dry mass term is shown in Figure 16. As expected, this term is of much lower significance 
than the added mass terms. 
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The load terms (other than drag) are added together and compared with the measured riser loads in Figure 17. The 
peak measured tensions are coincident in time with the peak added mass loads. The added mass peak magnitudes grow 
in significance with each stage, indicating the importance of modeling added mass for these canopies. As expected, 
the computed loads do not exceed the total riser tension at any point. Therefore, the difference between riser tension 
and these other loads must be due to drag area growth (neglecting spring elongation and apparent mass). 
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Figure 17. CDT-3-2 parachute load components compared with riser tension. 
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The isolated drag area for each canopy is plotted against enclosed volume in Figure 18. This information is used 
to determine key empirical look-up tables for each flight. The overall trend could be described as quadratic. Therefore, 
the data are fit for each stage with second order polynomials for simplicity. The first and second stages (inset) show 
concave-down trends. This is most likely due to the cluster interference effects at the end of the reefed stages. The 
large differences between curves for both canopies are most likely due to aerodynamic and surface shape differences. 
Full open data are only fit through the end times of the curve fit plots shown previously. Therefore, no attempt is made 
to match the full open data in the upper right which oscillates about a volume of 120,000 ft? and drag area of 10,000 ft. 
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Figure 18. | CDT-3-2 isolated drag area vs. enclosed volume. 
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The curve fits are then used to compute the drag area growth with time. The drag area is multiplied by interpolated 
dynamic pressure to compute the isolated canopy drag forces, plotted in Figure 19. Each stage shows a sudden drop- 
off in drag just after each disreef. This is most likely due to the sudden loss of pressure applied to the lower end of the 
canopy as the skirt suddenly expands. This sudden drop-off in load has been observed for some time, but has not been 
modeled in any CPAS simulation until now. 
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Figure 19. CDT-3-2 isolated drag force. 
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The curve fits are then used to compute the drag area growth with time. The results are plotted in Figure 20. For 
this case, the volume method succeeds in matching the observed peak riser tension timing and magnitude, especially 
in the disreef to full open. FAST simulations of disreef peak loads generally match the magnitude but occur earlier. 
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Figure 20. CDT-3-2 riser tension validation of peak load timing and magnitude (Kae = 0). 


In order to determine the sensitivity to apparent mass terms, the data reduction and simulation script was re-run 
for varying Ke values of 0.1, 0.25, 0.5, 1.0, and 2.0. Value of 0.1 and 0.25 (indicating that apparent mass is always 
10% or 25% of enclosed mass) had little effect on the resulting simulated riser tension. As Kae increased, the total 
added mass increased and the inferred drag area decreased, such that the sum of all terms was equal to measured riser 
tension. At about K,. = 0.5, the total estimated added mass for the lead canopy (Main S/N 6) exceeds the peak measured 
riser tension, making it impossible to properly solve for drag area. The lagging canopy (Main S/N 5) did not encounter 
such a discrepancy until Ka. exceeded about 1.4, indicating that apparent mass is dependent on canopy shape. 
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B. CDT-3-4 

CDT-3-4 was a 3-Main test where one of the canopies skipped second stage via a pre-cut reefing line.?° One of the 
non-skipping canopies was almost always obscured while the other was visible in the first stage and full open. The 3- 
D models are compared to the images in Figure 21. Due to the non-skipping canopies being crowded out, only data 
on the skipped-stage parachute, Main S/N 4, are available during the skipped stage. 


Figure 21. | CDT-3-4 selected 3-D canopy models of 1 stage and skipped 2™ stage. 
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The computed perimeter is plotted in Figure 22. The first stage perimeter for S/N 4 is consistently larger than the 
known reefing line length, perhaps indicating a bias in the scale factor used in the side view photo. Although the 
frequency of chase photos was limited, the data make a good transition to the KX photogrammetric results after each 
parachute disreefs to full open. 
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Figure 22. CDT-3-4 skirt perimeter. 
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The enclosed volume has a linear trend with projected area, as shown in Figure 23. This trend is consistent with 
all CPAS tests analyzed to date. Because the chase photo imagery is limited, this relationship allows for the use of KX 
projected area data to be transformed into volume data near the time of full open of each canopy. 
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Figure 23. CDT-3-4 enclosed volume vs. projected area linear trends. 
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The early chase photo-based projected area data were combined with later KX data to create composite projected 
area data time histories. These are plotted as pink and light blue circles in Figure 24. This also allows for using the 
more accurate KX data during later oscillatory periods. 
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Figure 24.  CDT-3-4 projected area. 
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The directly measured enclosed volume data are compared with the inferred composite data in Figure 25. The 
volume vs. time curves were fit with third order polynomials for all the stages of Main S/N 4 and S/N 6. 
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Figure 25. | CDT-3-4 enclosed volume. 
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The directly measured enclosed volume data are compared with the inferred composite data in Figure 26. 
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Figure 26. CDT-3-4 enclosed volume rate. 
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From there, analysis of each load term proceeded as in the previous example. The sum of each of the computed 
load components are compared with the measured riser tensions in Figure 27. As expected, there is a significant peak 
in added mass load coincident with the measured peak riser tension in the skipped second stage. The difference 
between the measured riser tension and computed load is assumed to be exclusively due to drag area growth. 
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Figure 27. CDT-3-4 parachute load components compared with riser tension. 
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The resulting drag area vs. volume conversion plot is presented in Figure 28. The data are again matched with 
quadratic curves. The curves for first stage (inset) are much smaller than the full open data. Oscillatory data after full 


open are not included in the matching process. 
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Figure 28. CDT-3-4 drag area vs. enclosed volume. 


The drag areas for both simulated canopies were then computed using the look-up table. Simulated riser tension 
are compared with the measured data in Figure 29. Again, the peak inflation timing and magnitude are an acceptable 


match. 
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Figure 29. CDT-3-4 riser tension validation of peak load timing and magnitude. 


As with the previous test, small increases to the amount of apparent mass had little effect on matching the measured 
riser loads. This was the case until Kae increased to about 0.6, at which point the simulated added mass matched the 
riser load, leading to the physically impossible situation of zero drag area. 


IV. Conclusion 


The current CPAS simulation tool (FAST) generally matches the peak load magnitude, but with early timing, 
implying a flaw in the added mass model. Examination of the current added mass algorithm shows many 
undocumented features which may be problematic in extreme conditions. 

A method was developed to describe the parachute inflation state using time-varying volume data, unlike 
traditional parachute simulation methods which calculate drag area growth based on empirical data. Volume data were 
obtained from high quality imagery of chase photos and upward-looking HD video of CPAS Main 116 ft Do ringsail 
canopies in clusters (described in a companion paper). The volume and volume rate data were then combined with 
best estimate atmospheric and vehicle state data to estimate parachute added mass components. When the sum of the 
load components are subtracted from measured riser tension, the parachute drag can be isolated. Drag area is modeled 
as a function of parachute volume. 

The curve fits were then combined with interpolated test data to simulate peak riser tensions, which match in both 
magnitude and timing. The reconstructions presented assumed that apparent mass due to viscosity was negligible, 
though additional data, such as CFD solutions for drag, could be used to better estimate this effect. 

By matching data to polynomial functions, this volume method may be difficult to incorporate into Monte Carlo 
simulations. Further, the curve fits are very sensitive to the data selection end point, which is subjective. This study is 
a preliminary effort to generate a fully documented added mass model relevant to large ringsail canopies anchored to 
high quality modern data. 
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